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We consider the hexagonal pattern forming in the cross-section of an optical beam produced by 
a Kerr cavity, and we study the quantum correlations characterizing this structure. By using argu- 
ments related to the symmetry broken by the pattern formation, we identify a complete scenario of 
six-mode entanglement. Five independent phase quadratures combinations, connecting the hexag- 
onal modes, are shown to exhibit sub-shot-noise fluctuations. By means of a non-linear quantum 
calculation technique, quantum correlations among the mode photon numbers are demonstrated 
" l". and calculated. 
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^ ' I. INTRODUCTION 

<\ 

, It is well-known that spatial patterns may arise spontaneously in the transverse cross-section of optical beams as a 
consequence of nonlinear wave mixing processes . Equally fascinating, by maybe less known, is the presence of highly 
' non classical spatial correlations in the beam cross-section, underlying the classical process of pattern formation. In 
^ optical systems the nonlinearity is associated with the simoultaneous absorption and emission of a number of photons, 
a circumstance that creates entanglement among the waves that form the pattern. On a macroscopic level, this is at 
the very origin of non-classical spatial correlations in the transverse far-field plane 0|b. 

Among the possible spatial structures that may spontaneously arise in a beam cross section, the hexagonal pattern 
is one of the most common ones. Hexagonal structures have been predicted to form, typically, in Kerr or Kerr-like 
T-H media in various configurations. These include counterpropagating waves interacting with a slice of Kerr material in a 
_ cavityless configuration Q , systems with a single feedback mirror j^] , and planar resonators filled with a self- focusing 
1 Kerr medium Q, or with a non-saturable absorber |^]. Experimental observations of this kind of optical pattern 
Q^l include counterpropagating beams in sodium and a variety of systems with single feedback (see e.g. 0). 
^ ' , , This paper is devoted to the analysis of the quantum properties of the hexagonal pattern, in the context of a model 
P5 ■ for a Kerr medium enclosed in a planar resonator driven by a plane-wave pump beam. This model represents, for 
^ , its simplicity, a paradigm for optical pattern formation. In fact, it was one of the the first models able to predict a 
tranverse modulational instability in an optical system Q . If only one dimension in the transverse plane is considered 
1^ , the modulational instability gives origin to a roll pattern in the near field immediately out of the cavity. The two 
dimensional version of the model, analysed by [Q, predicts that, above a suitable threshold of the input intensity, a 
spontaneous breaking of the translational symmetry in the transverse plane gives rise to a hexagonal pattern in the 
near field, which in the far field corresponds to six bright spots in hexagonal arrangement surrounding a central spot. 
Quantum aspects of the hexagonal pattern in the Kerr cavity were first evidenciated by Grynberg and Lugiato 
. Quantum correlations in the light intensities of groups of four among the six hexagonal spots were predicted by 
using momentum-energy conservation arguments, which did not take into account the effects of dissipation through 
the cavity mirror (hence valid only for a single pass through the nonlinear medium). 

As in the approach of we will restrict our analysis to a seven-mode model, valid close to the instability threshold. 
By using a quantum calculation technique which takes into account the full nonlinearity of the problem, we will be 
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able to demonstrate that the results of |10|. hold also in the presence of cavity dissipation. Moreover, we will derive 
an analytical formula for the fluctuation spectrum describing the quantum correlation among mode intensities. 

In the second part of the paper we will resort to more traditional calculation technique, namely to the usual 
linearization of the quantum fluctuation dynamics around the classical steady-state pattern. By using arguments 
related to the translational symmetry of the model, broken by the pattern formation, we shall identify a complete 
scenario of quantum correlation among the phase quadratures of the six hexagonal modes. 

We believe that the six mode entanglement described in this paper, besides being likely to be accessible to experi- 
ments, might be of relevance from the point of view of quantum information applications. 

Section || introduces the quantum mechanical and classical models describing the self focusing Kerr cavity. Section 
is devoted to photon number correlations among the hexagonal modes in the context of a nonlinear model. Section 
investigates in general the phase quadrature correlation among modes. Section |^ concludes and discusses open 
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questions. 
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II. THE MODEL 



In this section we present the quantum mechanical counterpart of a well-known semi-classical model , which 
describes the dynamics of the scalar electric field in a cavity filled with an isotropic Kerr medium. 

We consider a one-directional planar cavity (Fig. |l|), with a single input-output port, driven by a coherent, plane- 
wave, monochromatic field of frequency ujq. The model is derived in the framework of the slowly varying envelope and 
paraxial approximation, and of the cavity mean field limit |l2| , which allows to neglect the dependence of the field 
on the longitudinal coordinate z along the sample. Under these assumptions only one longitudinal cavity mode is 
relevant, precisely the one corresponding to the longitudinal cavity resonance ojc closest to ujq. We denote by A{x,t) 
the intracavity field envelope operator corresponding to this mode; it depends on the transverse space coordinate 
X = {x, y) and time t, and obey standard equal-time commutation relations 

[A{x,t),A\x',t)\^5{x-x'). (1) 

By adopting a picture where the fast oscillation at the carrier frequency ojq is eliminated, the reversible part of the 
intracavity field dynamics is governed by the three Hamiltonian terms 

H = Hfree + H^xt + Hint ■ (2) 

The coupling due to the Kerr nonlinearity in the medium is accounted for by 

Hrnt - -?i7| j d'x [A\x)f [A {xf (3) 

where the coupling constant g is proportional to the third order x^^'^ susceptibility, and 7 is the cavity linewidth. This 
Hamiltonian describes four-wave-mixing interactions at a microscopic level, where they correspond to simultaneous 
annihilation and creation of photons in pairs Free propagation in a planar cavity is described, in the paraxial 
approximation, by 

Hfree = f d^xA\x) (A - Zf^V^) A{x) , (4) 



where A is the cavity detuning parameter; the two dimensional transverse Laplacian = + models the effect 

of diffraction in the paraxial approximation; the parameter Ijj « a/ LA/(27rT), with A being the light wavelength, L 
the total cavity length, and T the transmittivity coefficient of the cavity input-output mirror, defines the characteristic 
length scale for transverse pattern formation. Finally, the coherent pumping, by a classical plane-wave driving field 
of amplitude E'in , is modeled by the term: 

Hext = «?i7 / d'^x [£;,„At(f) - E*nA{x)] . (5) 



Dissipation through the cavity mirror can be described in the framework of standard cavity input-output formalism 
(see e.g. jl^). As a result the dynamics of the intracavity envelope operator is governed by the equation 

^A(f , t) = -7 [(1 + zA - iiys/^) A{x, t) - igA{x, t)^A^{x, t) - E,„] + ^SAn^x, t) , (6) 

where the noise operator SA{x,t)in represents the vacuum fluctuations entering the cavity input-output mirror. It 
has zero mean, and obeys the free field commutation relation 

M,„(f,i) , ML(f' ,t')] = S{x-x')5{t - t') . (7) 

Equation (^ must be coupled with the relation at the cavity mirror, linking the outgoing reflected fleld Aout{xjt) 
with the intracavity and input fields 



Aout {x, t) V 27^(^, i) - A^n {x, t) , (8) 

where, in our notation, the total input field operator is given by Ain{x, t) = ^^l2Ein + bAinix, t). 
The classical counterpart of this model is the well known partial differential equation |9| 
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1 OEJ 

-— = -E + E,„ + i{\E\''-A + llV^E)E, (9) 

for the envelope of the classical electric field E. This equation is straightforwardly obtained by its operatorial version 
(O) by dropping the noise term, and scaling the input and intracavity field intensities to the saturation photon number 
f79,i.e. g\E\^^\E\^,g\E,^\^^\E,^\\ ^ ^ 

Since the model has an overall translational symmetry in the transverse plane, equation admits a transversely 
homogeneous stationary solution Eqs, which obeys the well-known cubic steady-state equation |14[| 



'in\ — \-t^Os\ 



l + {A~\Eos\^f . (fO) 



As shown by Lugiato and Lefever by increasing the input intensity a spontaneous breaking of the translational 
symmetry occurs, and the steady state (^o|) becomes unstable with respect to perturbations spatially modulated in 
the transverse plane. As it is common in pattern formation processes ||l| the instability develops with a definite critical 
wave number fee, which characterizes the periodicity of the pattern immediately above the instability threshold. The 
critical point for the onset of the modulational instability is given by 

\Eos\^ = l, llkl = 2-A. (11) 

Among the transverse modes on the critical circle |fc| = kc, the Kerr nonlinearity selects a discrete set of six wave 
vectors |^,|l^ , so that immediately above the critical point the near field intensity distribution consists of a hexagonal 
lattice of bright spots , while the far field shows a central bright spot surrounded by six spots in hexagonal arrangement. 

The origin of this particular pattern can be understood by considering the microscopic processes leading to off-axis 
emission of photons allowed by energy-momentum conservation [|6| . By referring to the scheme of figure a first 
process that preserves the radiation transverse momentum involves destruction of two photons of the plane- wave input 
beam, and creation of two photons propagating off-axis in symmetrical directions. This is the only possible process in 
a 1-D version of this model and would lead to the formation of a stripe pattern in the near field, as e.g. predicted 
by a vectorial model for a self-defocusing Kerr medium [OJlSl, or in the case of a degenerate optical parametric 



oscillator 19 1. In this system, however, the stripe pattern is unstable, because secondary four- wave mixing processes 
can take place. By referring to figure ^, we can e.g. consider destruction of one on-axis photon of the input beam and 
one photon propagating in direction 1, and creation of two off-axis photons propagating in directions 2 and 6. By 
considering also the symmetric process, that involves destruction of a photon in mode 4 and one in the on-axis mode, 
and creation of two photons in mode 3 and 6, respectively, we see that the six off-axis hexagonal modes are activated. 
Close to the threshold for pattern formation, all the other processes allowed by energy-momentum conservation, that 
do not involve on-axis photons, can be neglected, because off-axis mode intensities are much weaker than that of the 
on-axis homogeneous mode, where the pump beam is injected. 



III. PHOTON NUMBER CORRELATION IN THE HEXAGONAL PATTERN 



The microscopic processes which generate the hexagonal pattern are at the origin of non-classical correlation among 
different portions of the far-field cross section. Referring to the example of figure ^, it can be noticed that whenever 
two photons are detected at positions 2 and 3 in the far field, there must be two other photons at position 5 and 6, 
which suggests the existence of a high level of correlation between the sum of photon numbers N2 + A3 and + Nq 
(the same reasoning can be applied to any symmetric group of four spots). 

A more quantitative argument for the existence of such a kind of intensity correlation in the hexagonal pattern was 
given by Grynberg and Lugiato [l0|] . They considered a discrete model for the hexagonal pattern formation in Kerr 
media, which takes into account only the six hexagonal modes plus the homogeneous mode = 0, valid close to the 
threshold for pattern formation, where these modes gives the most relevant contribution to the system dynamics. By 
considering a quantization box in the transverse plane of side & (6 ^ 00 at the end of the calculations) they set 

1 1 

A{x,t) = -ao{t) + -Y^a,{t)e^'^-^ , (12) 

with (see figure |b) fci = (0, /cc), k2 = (fcc/2, V3/2fcc), h = (fcc/2, -\/3/2fcc), k4 = -fci, fcs = -k2, ke = -k^. By 
inserting this expansion into the system Hamiltonian (|^) they found that the combination of mode photon numbers 

= AT, + A,ei - (A^^ffi3 + A,e4) , i = l,...,6, ^^©j = ^+Jmod6 (13) 
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with Ni = a|ai, commutes with the discretized Hamiltonian. Thus, in a single pass through the nonhnear medium 
the observable N- is preserved. If there are no photons in the hexagonal modes at the crystal entrance, perfect 
correlation is found at the exit between e.g. Ni + Nif^i and Niq)^ + Ni(^4. However, when the crystal is inserted in 
an open resonator, the field is recycled by the cavity and dissipation through the cavity input/output mirror should 
be taken into account. As a matter of fact, the release of a photon out of the cavity mirror is a random process 
uncorrelated with the nonlinear processes taking place in the crystal, so that one would expect that cavity dissipation 
is detrimental for this kind of correlation. 

However, we will now show that this is not the case. In order to elucidate this point, we shall follow two different 
techniques. The standard one, which will be employed in the next section, makes use of the small quantum noise 
approximation, and consists of linearizing the dynamics of quantum fluctuations around a classical steady state. The 
other one consists of trying to solve the non- linear Langevin equations that govern the dynamics of the photon number 
operators of the hexagonal modes. These equations take the general form (see e.g. 



dt 



-N,, 



1 



(14) 



with a]"" being noise operators modelling the vacuum input fluctuations on i-th mode. They obey standard free-field 
commutation rules 



and have the noticeable commutator with a generic intracavity operator O jT^ 



(15) 



(16) 



where 9 is the step function. Equation (|16|) expresses the fact that intracavity operators at time t commute with the 
input operator at a later time t' due to causality. The discretized version of Hamiltonian (^, obtained by means of 
the expansion (|l|), results 



r(i) 



(2) 



(3) 



Hint — HsPM + -ffcPM + ^^FWM ^ ^^FWM ^^FWM ' 

where i?spM describes self-phase modulation 



3=0 



while HcpM describes cross-phase modulation 

HcPM = -2fi-fgy^alaia]a 



' 



(17) 



(18) 



(19) 



and i?FWM describes four-wave mixing 



(1) 

FWM 



2<J 



h.c. 



(20) 



(2) 

FWM 



^ a^ai©3a]a]g)3 + h.c. 



(21) 



(3) 



'FWM - -2?i75 

Other terms of the Hamiltonian commute with the photon number operators 



(22) 
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The second term at RHS of the dynamical equations (|Tj) is therefore a comphcated expression, involving third 
order products of mode operators, so that finding an explicit solution seems unlikely. However, when considering the 
combination of mode photon numbers -/V_ given by Eq.(O) the dynamics drastically simplify: 



j^N-=^2^N--'-[N-,H] + ^G-, (23) 
= -277V_ + , (24) 



where we used the fact that N_ commutes with the system Hamiltonian, and the operator is defined as: 

~ ("les^ies + alesOiea) - (ale4"ie4 + '^Im'^ieij (25) 

When the input operators of the hexagonal modes are in the vacuum state, (G_(t)) = 0. Moreover, by using 

i) the commutation relations of input operators, 

ii) the causality relations expressed by Eq. (pj|), 

iii) the fact that the a™ annihilates the vacuum on the right and a| ™ the vacuum on the left, 
it is not difficult to show that 

{G.{t)G.{t'))^S{t-t'){N+{t)) , (26) 

with 

N+ = N,+ iV.ei + ^^©3 + N,(S4 . (27) 
At steady state (iV+(t)) does not depend on time and Eq.(p4|) is easily solved in the frequency domain. By setting 

iV_(c.)= I -^7V_(i)e-', (28) 

we have 

7V_(w) = V2^-^^^. (29) 

By using the boundary relations for the outgoing mode operators a""* = — y/2r/a™', we get 

7V°"* (w) = 27iV_ (w) + TV™ (w) - y^G- (w) = V2^G_ (w) „ + TV!" (w) . (30) 

27 — iuj 

The fluctuation spectrum of the output photon number combination N°_"* can now be calculated as: 

V^u;) ^Jdt (^7V°"*(t)iV™*(0))e-* = 7(^+> - ^J^^) , (31) 

where 5N'^^ — N°^'-* ~ {N'i"*) — iV"""* . Hence the noise spectrum of A^°"* turns out to have a standard Lorentzian shape 
[ pT| , with the Lorentzian dip well below the shot-noise level represented by 2^(N^) = (A^™*) . This result is reminiscent 

of the ' twin beam" correlation spectrum calculated for the down-converted fields generated by a parametric oscillator 
without spatial aspects pO| . As a matter of fact, in the case of twin beam generation, the photon number difference 
between the two beams has zero commutator with the system Hamiltonian. The calculation of the noise spectrum of 
the number difference may be performed following exactly the same steps outlined here, and it is not a coincidence 
that the result turns out to be the same. However, in contrast to the more standard method of analysis, which exploits 
the small quantum noise limit by linearizing the model equations around a classical steady-state, the method used 
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here is "exact"[J, because it takes into account the full nonlinearity, and does not depend on the particular steady 
state chosen. 

In this connection it should be noted that our result concerning spatial intensity correlation in the hexagonal 
pattern is valid both above and below the threshold for pattern formation, where there is no classical pattern at all. 
Off-axis emission of photons is here generated purely by quantum fluctuations; on average translational and rotational 
symmetry are preserved, so that off-axis photons can be emitted in any direction. However, the mechanism of photon 
absorption and emission is such to preserve the transverse radiation momentum; as a consequence any six regions 
in hexagonal arrangement with arbitrary orientation in the far field plane display exactly the same kind of intensity 
correlation that we have described above threshold. 

Incidentally, it is worth noting that the same method can be applied to analyze the quantum spatial correlation 
in the far field plane of the down-converted field emitted by a optical parametric oscillator, where under proper 
circumstances the downconverted field is emitted in the form of two plane waves slightly tilted with respect to the 
cavity longitudinal axis [ p9| . By interference these two waves give rise in the near field to a stripe pattern, while the 
far-field intensity distribution consists of two symmetrical spots. From a microscopic point of view, the most relevant 
three-wave-mixing process close to threshold involves destruction of a pump photon that propagates longitudinally, 
and creation of two down-converted photons propagating in symmetrical directions, as it is required by transverse 
momentum conservation. Hence, the two off-axis signal waves emitted above threshold are made of twin photons and 
the photon numbers iVi , A'2 that cross in the unit time two ideal detectors surroundings the two spots in the far field 
plane are highly correlated pl|-|2^ . In the framework of a three-mode-model, valid close to threshold, a calculation 
analogous to that performed in the case of hexagons shows that the noise spectrum of their difference ~ N2 — Ni 
is given exactly by the formula 

A last interesting remark is the following: in the Kerr cavity model, when the quantum fluctuation dynamics 
is linearized around the homogeneous steady-state below threshold, intensity correlation is found between any two 
symmetric wave- vectors k and —k close to the critical circle, just as in the case of parametric down-conversion. Hence 
the linear analysis does not give any hint about the existence of a hexagonal pattern above threshold, with its quantum 
correlations. This is not in contradiction with our result, since quantum correlation between two symmetric modes 
obviously implies a noise reduction in the observable 7V_ connecting two couples of symmetric modes. However, the 
inverse is not in general true. Since the intensity difference between symmetric modes is not a constant of motion of 
the full Hamiltonian, we expect that the level of correlation between symmetric modes decreases when approaching 
the region where linearization fails, that is, in the neighbourhood of the instability point, where the size of fluctuations 
increases, or for a truly microscopic system characterized by a small saturation photon number parameter. On the 
contrary, our non-linear analysis shows that the correlations which are at the origin of noise suppression in iV_ do 
not depend on the distance from theshold or on the system size. 

IV. PHASE QUADRATURE CORRELATION IN THE HEXAGONAL PATTERN 

In the previous section it has been demonstrated a high level of intensity correlation in groups of four among the 
hexagonal modes. A natural question which arises is whether there exist any other kind of correlation among the 
six modes at the quantum level. This possibility is suggested by the example of twin beams generation. In fact, in 
this case, not only the photon number fluctuations are correlated, but also the "phase" fluctuations of the two beams 
are anti-correlated at a quantum level. To be more precise, there exist two orthogonal field quadratures of the two 
beams which show, at the same time, a high level of (anti)correlation. This is at the origin of the EPR aspects of the 
twin beams p5|j2^ , which are widely exploited in the field of quantum information with continuous variables (see e.g. 

&!§)• 

Since we do not have in mind other constant of motion (other observable commuting with the Hamiltonian) jWe 
have to resort to traditional means of calculations. Namely, in the framework of the seven-mode discrete model 
we separate the quantum fluctuation operator from the classical mean field 

ai{t) ^ ^/3, + Sttiit) i = 0,l,...6 (32) 



'^More precisely, it is exact as long as transverse modes other than those forming the pattern at the instability onset are 
negligible. 
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where (3i are the classical steady state amplitudes of the seven mode. More precisely, they are the steady state 
solution of a set of classical dynamical equation, obtained by introducing in Eq.(|p) the discrete expansion of the 



envelope operator E{x, t) 
1 d 



jao(t) + \ oti{t)e'^^^'^ . These equations have the form: 



1 d 
^'dt 



ao^Ei- (1 



/ J 



= —(1 + 2i)aj + ioLjO^j + 2jQ;j|aop + 2iQ; 



6 



,6, 



(33) 



where we took into account the fact that all the hexagonal modes have the same critical transverse wave number, 
such that A + ljjkl ^2. 

By inserting the expansion into the model equation (^, with the ansatz (^2[), and keeping only the leading 
terms in the small quantum fluctuations we are left with a problem of the form 



d_ 

dt 



foe 



I 



14 X 14 
matrix 



matrix elements 
are functions of 
V(/3o,/3o*...,/?6,/36*)/ 



fog 
foj, 



for 

fo; 



Vfot" 



(34) 



where the explicit form of the matrix elements in terms of the classical steady state amplitudes (/3o, /Sq • ■ ■ , /?6j /Sg ) is 
given in Appendix A. The input fluctuation operators at RHS are in the vacuum state and have commutation relations 
as those in Eq.jisl). This represents a 14 x 14 linear problem, which is trivial from a numeric point of view; however, 
when searching for some explicit combination of mode operators that has sub-shot-noise fluctuations, the problem 
seems too complex to find analytical solutions. Nevertheless, some hints came from the analysis of the classical steady 
state. 

We integrated numerically the set of classical equations (^) , and looked at the long time behaviour of the mode 
amplitudes. Above the critical point (O) the hexagonal steady-state of the classical equations of the model 



1 1 ^ - 

is characterized by: 

- All the hexagonal modes have the same mean intensity 

|/3i| = 1/321 = ... 1/361 := |/?| 



(35) 



(36) 



so that we set (ij — e"^^ . Figure |^a shows the steady-state hexagonal mode amplitude as a function of the 



input field intensity. The figure is obtained by integrating the classical model equations (|3^) under a slow (with 
respect to the characteristic time needed to reach the steady state) increasing of the input field intensity across 
the instability threshold (solid curve); the dashed curve correspond to a slow decrease of the input intensity. 
As it is well know [p|, p^ the instability is subcritical, and the hexagonal mode amplitude shows the typical 
hysteresis cycle. Figure^ is the same for the amplitude |/3o| of the homogeneous mode. 



The sums of the phases of symmetric modes are all equal 



(37) 



The differences of the phases of symmetric modes sum up to zerof 



^Notice that in the classical analysis presented in some initial assumptions on the steady state were made so that 

condition (|37|) was automatically fulfilled,with = 0, and condition (Bq) was written as (Ai^i 4- A^a -I- A<^5)/2 — 0, tt. 
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A(/3i + A</?3 + A<p5 = with Atfj = tfj — (^j^3 . 



(38) 



There are two phases which are not fixed by the steady-state equations, namely two of the phase differences between 
symmetric modes, say and A933. The value of these phases at steady-state depends only on initial conditions, and 
these variables are dominated by noise. This circumstance is a noteworthy consequence of the translational symmetry 
of the model, broken by the formation of the pattern in the transverse plane. As a matter of fact, for each value of 
the input field intensity \Ei\'^ above the critical point, there exists a continuous set of steady-state solutions of the 
model equations (p3), of the form 



icj) 



COS (^ki • X + A(pi^ + cos (^ks • x + Aip^^ + cos • x — {A(pi + A(^3 



(39) 



corresponding to an arbitrary choice of A(/?i, Ai^3. By setting Aipi = ki ■ Ax, and A(p3 = ^3 • Aa:, we notice that 
Aipi + A(f3 = (fci + fc3) • Ax — — ^5 • Ax. Thus the set of solutions (|3^ ) corresponds to the continuous set of rigid 
translations of the hexagonal pattern in the {x, y) plane, x — > x + Ax. 

As it is common in continuos symmetry breaking phase transition pO| ] close to the critical point, the noise is con- 
centrated on the mode that aims to restore the symmetry broken by the transition. In our case we argue that close 
to the instability threshold, quantum fiuctuations corresponding to rigid translations of the pattern result in huge 
fiuctuations of the differences of the phases between symmetric modes A^pj , but leaves invariant the quantities 



{A) Aifii + Aifis + Aip5 

{B) (.^, -I- .^,;e3) - (<Pi©i + <Pj©4) i = l,...6 

Hence, our possible candidates for squeezing are the observables 

{5ai — 5a4 -I- Sa^ — Jag + 5a<^ — 5a2) e~'^ -I- h.c. 

which corresponds to the classical quantity (A) in ( pO[ ) for 6 = (j) + n/2, and 

[{Sa^ + 6ai(S3) - {Sa^Qi + Sai^^)] e""'' + h.c. , 



(40) 
(41) 



(42) 



(43) 



which corresponds to (B) in ( ^l| ) when t/j = 
optimization parameters. 

In addition, an obvious observable to take into account is 



7r/2. In these definitions 0,ip are left arbitrary and will be used as 



[{Sai - (5aj©3) + (fci©i - StttQ^)] e"^ 



i = l,2, 



(44) 



In the small quantum noise approximation, this observable corresponds to the the photon number difference A^_, 
which was shown in the previous section to have a sub-shot noise fiuctuation spectrum. 

Luckily enough, when considering the combinations of mode quadratures given by ( ^ ) , ( 43 ) , ( ^ ) , their dynamical 
equations decouple from the other equations of the system (|3j), and we are left with 2x2 linear problems. 

Before examining specific cases, we note that the quantum dynamics of a generic linear system can be written as 

d 



dt 



V = MV 



27V" 



(45) 



where V, V" are the system and noise operators vectors respectively, while M is a coefficient matrix. 

Suppose that the components of the vector V are two conjugate field quadratures Z{Q) and Z{ti/T}\ at steady state 
the system can be solved in the frequency domain, and, with the aid of an input-output relation 

(46) 



it is possible to calculate the output correlation matrix. In the frequency domain it reads 

C™*(w) = / da''(V°"*(w)[V™*]^(w')> = \2-i{m + iujiy^ +1\ C'"(tj) [27(M-iu;I)"^ -hi 



(47) 



where I is the identity matrix, T means the transpose, and C'"(a;) — J duj'{V"^{uj) [V'"]'^(a;')) is the input correlation 
matrix. 

Moreover, the noise spectrum of a generic quadrature 
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= Z{Q) cos^j + Z{Tr/2) sini/', 

is given by 

Sz{i^,u;) := J da;'(,5Z™*(V,w)<5Z™'(^,t^')> 

= C°"*(tj) cos^ V + C°"2*(w) sin^ V + (C°;'2*('^) + C™! (w)) sinV' cos V . 
With the above in mind wc arc now going to consider specific cases. 



(48) 

(49) 
(50) 



A. Noise in the sum of phase differences 



Let us consider the quadrature operator 



1 



W{9) = —j= (ai — 04 + 03 — oe + as — 02) e ' + h.c. , 
The dynamics of the fiuctuation vector V = [(5W^(0), (5W^(7r/2)]-^, gives rise to a closed system like (p^), with 



M = 7 



Re{^+} -Ini{y^+} 
Im{y4_} Re{yl_} 



(51) 



(52) 



and 



A± = [-1 -2i + 2i|/3oP + 10i|/3p - 4ipof3* - 4i/3o/3] 



(53) 



In writing ( |53| ) the symmetries of the steady state have been taken into account; moreover, among the possible 
stationary states (p9) we used the one corresponding to Aipi — 0, A(^3 = 0. The input correlation matrix is given by 



C" = 



1 i 

-i 1 



(54) 



We performed calculations of the noise spectrum Sw{(^,^) for various values of the input beam intensity |i?mp in 
the region where the hexagonal solution exists. The amplitudes of the homogeneous mode /3o, and of the hexagonal 
mode (3 were obtained both by numerically integrating the dynamical classical equations (33) and looking at the 



long time behavior, and by numerically solving the nonlinear steady-state equations (see Appendix B). An excellent 
agreement between the two approaches was found. 

Figure ^ shows an example of the typical results. Part Al of the figure plots the zero frequency spectrum as a 
function of the quadrature angle 9, shifted by tt/2. Sub shot-noise fluctuations are present for a rather broad range 
of quadrature angles; it has to be noticed that the quadrature operator exactly corresponding to the classical sum 
of phase differences (40) is the one with 9 + 7t/2 = cj) where is the steady-state phase of (5 (indicated by arrows 
in the figure). Hence sub shot noise fiuctuations are present for quadratures somehow rotated with respect to the 
semiclassical "phase" quadrature. In Part A2 of the figure the quadrature angle is chosen to optimize squeezing; the 
plots evidenciate the typical Lorentzian shape of the fiuctuation spectrum. 

Line (A) in figure plots the best squeezing (that is, the low frequency noise with the phase 9 optimized) as a 
function of the input field intensity, and shows that sub shot noise fiuctuations for this observable are present in the 
whole region where the hexagonal solution exists. 



B. Noise in the difference of phase sums 



We consider now the quadrature operators 



1 



Q(^) = (oi 

V4 



-I- aj©3 - aj©i - ai©4) e + h.c. , 



(55) 



In fact there are two independent mode combinations of this form, e.g. for z = 1, 2, and for both of them we get a 
closed system like (45) for the fluctuation vector V — [(5(5(0), (5Q(7r/2)]^, with 
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M 



Re{B+} 
lm{B-} 



-lm{B+} 
Re{B-} 



(56) 



and 



B± = [-l-2i + 2i|/3o|^ + 6i|/3|2- 



(57) 



The input noise correlations are again given by (^). Fig|| shows the noise spectrum SQ{ip,ijj), and is analogous 
to Fig.^ As in the previous case, the spectrum has a Lorentzian shape and exhibits squeezing in the whole region 
where hexagons are predicted by the seven mode model, although the squeezed quadratures do not exactly coincide 
with the one corresponding to the classical phases in (p|). 

Line (B) in figure ^ plots the best squeezing as a function of the input field intensity, and shows that also in this 
case sub shot noise fluctuations for this observable are present in the whole region where the hexagonal solution exists. 



C. Noise in the sum of intensity differences 



Finally, let us consider the two independent quadratures 



1 



X(V') = -y= [(fli + flies) - (ai©i + an^i)] e + h.c. i = 1, 2 



Again, their fluctuation vectors V = [SX (0) , 6X (t: /2)]'^ , give rise to two closed systems like ( p5| ) with 

M = 7 



Rc{C+} -lm{C+} 
Im{C-} Re{C-} 



(58) 



(59) 



and 



C± 



= [-1 - 2i 
± [2z/3*/3* ^ 



-2i\po\^ 



2i(3of3* + 2i(3*f3] 



(60) 



The input correlation matrix is the same of Eq. (^4|) . 

It is worth noting that the quadrature corresponding to the hexagonal mode stationary phase <p satisfies a simple 
equation 



(61) 



-<5X(0) = -2^5Xi^) + V27<5X" 



which gives a Lorentzian shape spectrum 



Sx{(f',uj) 



4^2 



(62) 



Taking into account that now the shot-noise level is scaled to 1, this is in perfect agreement with the results of Section 
II, since 5X{(f)) represents the linearized version of the observable N_. 

It is interesting to remark that, differently with respect to the previous cases, for this combination of modes the low 
frequency noise reduction abruptly disappears when the quadrature is slightly rotated with respect to the amplitude 
quadrature. In fact, as shown by figure ^, a large positive peak located at zero frequency appears on the Lorentzian 
spectrum when the quadrature angle is shifted from <f>. 

Line (C) in figure |^ plots a numerical evaluation of the best squeezing as a function of the input field intensity. As 
it is obvious from the analytical formula |6^ ) complete noise suppression is predicted at zero frequency for — cj) (the 
amplitude quadrature), in the whole region where the hexagonal solution exists. 
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V. CONCLUDING REMARKS AND OPEN QUESTIONS 



In conclusion, we have studied the quantum features underlying the hexagonal pattern formation in a Kerr cavity. 
We have identified a rich scenario of purely quantum correlations among the off axis modes that form the hexagonal 
pattern. In fact we have shown that at least five independent combinations of mode observables exhibit sub shot-noise 
fluctuations. 

We believe that the six- mode entanglement analysed in this paper might be of great interest. On the one side this 
kind of patterns are widely studied in labs, so that it should be easy to look for the described quantum effects. On the 
other side, the quantum properties of such systems make them good candidates for quantum information processing 
with continuous variable p7|- p9| , pT[ . In particular, the existence of several observables showing correlations at the 
quantum level opens the possibility of distilling a considerable amount of entanglement in practical situations, and to 
use it. 

In this connection an interesting question is whether this kind of spatial entanglement is characteristic of the 
particular model analysed here, or it is in general typical of the hexagonal pattern, independently of the detailed 
mechanism that underlie its formation. Further investigations about models closer to experimental situations where 
hexagons heve been found are hence in order. 

An other interesting point concerns the possibility that some squeezing is missing from our picture. The symmetry 
of the problem suggests the possibility that quantum correlations - more probably anti-correlations-exist between the 
the sum of all hexagonal modes and the pump mode. This would be closely similar to the situation analysed in p2| , p3| , 
in the context of a self-defocusing Kerr cavity, where quantum anticorrelations were found between the pump on-axis 
intensity and the intensities of two symmetrical off-axis modes. 

A final open question is whether these quantum correlations survive when including in the model the full continuum 
of transverse modes. As a matter of fact, being the bifurcation subcritical, the hexagonal mode emerges at threshold 
with a finite amplitude. Hence higher order spatial harmonics might have a non-negligible amplitude evening in the 
neighbourhood of the threshold, and photon scattering in these additional modes might spoil our correlations. This 
will be the subject of future numerical investigations. 
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APPENDIX A 



In this Appendix we give the explicit form of the linear system (^4|) for the dynamics of fluctuations of the hexagonal 
modes. In writing these equations we took into account the symmetries of the hexagonal steady state expressed by 
conditions (^),(p7|) and (^); moreover among the possible steady states in ( |39| ) we choose for simplicity the one with 
A(^j = , j = 1, 3, 5 (this correspond to fixing the origin of the coordinate axes in the transverse plane). We get 



(63) 



^^Saa = [-1 -ie + 2i|/3oP + 12z|/?P] Sao + [iPo + 6*/?^] ^4 

6 ^ FT 

+ [2i(3Ql3* + 2if3*(3 + 4i\f3\^]Y^ Sa^ + [22/3o/3] H + V ~ '^"o" ' 

-j^Saj = [2i/3o/3* -t- 2i(3l(3 + 4i|/3|2] 5aQ + [2i/3o/3 + 2i(3^] Sal 

+ [-1 - 2^ -f 2^|/3oP + 12i|/3p] Sa^ + [ifi^] 5a] 

+ [2iM* + 2iP*o(3 + Ai\l3\^] Saj^i + [2tl3^] Sa]^, 

+ [4z|/3|2] Saj^2 + [2tPof3 + 2tp^] Sa]^^ 

+ [2i\P\^] Sajes + [tPl + 6»/32] 

+ [4i|/3|2] 5a,e4 + [2z/3o/? + 2iP^] Sa]^^ 

+ [2i0oP* + 2iP*f3 + 4z|/3|2] Saj^5 + [2*/?'] + ' i = 1, • ■ • , 6 . (64) 
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The elements of the matrices M can be built up by appropriately combining some of the above equations. 



APPENDIX B 



This appendix gives the explicit expression of the nonlinear equations for the classical hexagonal steady-state. These 
are obtained by equating to zero the l.h.s of Eqs.(^); in order to find a numerical solution wc transformed them in 
four nonlinear equations for the real variables {uq, vo,ui,vi} introduced as 

I3q = Eqs + Eqs {uq + ivo) , (65) 
(3^^iu,+^v,) . (66) 

The resulting equations are 

= -Uq + (A - \Eas\'^)va - |£^osP ^uqVq + uivi + MqUq + ^VquI 

1 3 1 \ 

zujvi + UQUivi +v^ + -vovf + — 7=Wi , (67) 



2^3 ' " ' " ' 2""-'^ ' 2V3 ^ 

( 3 1 

= -t;o - (A - \Eos\^)ufi + |£;os|^ 2uo + 3u^ + -ul + vl + -vl + ul 

+ + "o^^o + ^"ow? + vqUiVi + -^^j=uivl^ , (68) 

/ 2 1 

= — Ul + 2vi — \Eos\^ vi2uqVi + 2i'oMi H — i=uivi H — j=voul 
y V3 v3 

2 5 2 2 2 / — 2 ^ 3 I / \ 

+ 2uoVoUi + UqVi + -^u^vi + —^uquivi + SvqVi + vSvqVi + -v^ j , (69) 



= 



-Ul — 2ui -I- |i?osP ^3mi -I- GuqUi + VSuf + 2voVi + "^^i 



5 1 2 5 \ 

-I- 3uqUi -I- V^uquI + -uf + 2uqVqVi H — pUo«i + VqUi H — —vqUiVi + -rUivf . (70) 
4 v3 v3 4 I 

The above equations are solved numerically with the help of Mathematica, by using, for simplicity, the condition 
A = |£'osP, which assures that Ein = Eqs (see Eq. dlQ). The solutions are compared to the steady-state mode 
amplitudes obtained by numerically integrating the dynamical equations (p3|); a perfect agreement between the two 
approaches has been found. 
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VI. FIGURES 
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FIG. 1. Scheme of a ring cavity containing a Kerr medium (K). The longitudinal and transverse axis are shown for both 
input and output field. 
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(a) 



Fig. 2 



(b) 



FIG. 2. (a)Scheme of the microscopic processes leading to hexagonal pattern generation (see text), (b) Transverse wave 
vectors forming the hexagonal pattern. 
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FIG. 3. Example of the steady-state curves for the hexagonal mode (a) and the homogeneous mode(b) amplitudes as a 
function of the input beam intensity. Solid (dashed) lines are obtained by a numerical integration of Eqs.(^3|) over times long 
compared to time scales of dynamical evolution, and by a slow forward (backward) sweep of the input beam intensity across 
the instability threshold — 1). Cavity detuning is chosen as A = |i5osp. 
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Fig.4 



FIG. 4. Part (Al) shows the zero frequency spectrum of the observable W{6), given by Eq.(|5l[), vs the quadrature angle 
8 shifted by 7r/2, for two values of the input field. The arrows indicate the corresponding hexagonal mode phases. Part (A2) 
shows the frequency spectrum for the optimal value of quadrature angle. A — l-BosP- Other parameters are indicated in the 
figure. Dots correspond to a Lorentzian fit. 



17 




hex. mode 
phase 



(B1) 



IE_j'=0.97 

IE_j'=l.l 



1,5 2,0 2,5 3,0 

\|/ - *a (rad) 




Fig.5 



FIG. 5. Part (Bl) shows the zero frequency spectrum of the observable given by Eq.(^5|), vs the quadrature angle 

-0 shifted by 7r/2, for two values of the input field. The arrows indicate the corresponding hexagonal mode phases. Part (B2) 
shows the frequency spectrum for the optimal value of quadrature angle. A = |iSosp. Other parameters are indicated in the 
figure. Dots correspond to a Lorentzian fit. 
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FIG. 6. Noise spectrum of the observable X{ip) given by Eq.(^8|). The solid line is the Lorentzian spectrum when the 
quadrature angle ip = (f); the other lines correspond to quadratures slightly rotated with respect to the amplitude quadrature, 
as indicated in the figure. 
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FIG. 7. Best squeezing (low frequency noise with t he quadr ature angle optimi zed) a s a function of the input field intensity, 
for the three observables considered in Sections [V A (line A), IV B (line B), and IV C (line C). 
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